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ABSTRACT 


A problen of routing earth resource survey 
aircraft, proposed by NASA, is formulated as a 
traveling salesman problem, in which the salesman 
(aircraft) has a range constraint. A heuristic 
algorithm is presented, which seeks a minimal length 
set of subtours through n cities. The aircraft begins 
a subtour at the base location, visits a subset of the 
nm cities and returns when the range constraint 
prevents a visit to another city. Additional subtours 


are created until all cities are visited. 


The algorithm 1s programmed in FORTRAN for use on 
digital computers. The IBM 360/67 computer at the 
Naval Postgraduate School was used to find solutions 
to three operational problems of size seven, eighteen 
and twenty-five cities. Computation times for each 
problem was under 20 seconds and the solutions were 
Significiantly better than feasible solutions 


calculated without the use of the algorithn. 
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The National Aeronautics and Space Administration is 
fostering many technology utilization programs to explore 
the application of its aerospace technology to other fields. 
Suggesting how our finite natural resources can be used most 
effectively is one of itS more important considerations. 
NASA has been working to sophisticate the remote sensing of 
the earth from a variety of aerospace platforms. The 
information available from such a program can provide 
accurate and timely data to assist in the management of the 
earth's natural and cultivated resources. NASA has 
demonstrated the usefulness of remote sensing to a wide 
variety of users with emphasis on the disciplines of 
geology, agriculture, forestry, oceanology, meteorology and 
hydrology. For example, an agriculturist can be provided 
With identification and area measurements of major crop 
types for use ln forecasting potential crop yield. In 
addition, remote sensing can sometimes detect agricultural 
diseases and pest attacks before they can be detected fron 
the ground. More detailed discussions of the type of 
instruments used and examples of the data available from 
remote sensing can be found in references 4, 6, 7, 13, 15, 
and 16. 


In recent years land has emerged as the one resource 
"basic to both the good life and the goods of life" [ref 8}. 
There has been emphasis on efficient land use planning, 
particularly in metropolitan areas. Many experts are 
convinced that a well conceived program of remote sensing of 
metropolitan areas will enhance the quality of urban design. 


This program could provide, for a given area of land, its 





biological productivity, soil limitations, topology, access 
to transportation and utilities, existing vegetation and the 
use of adjacent land areas. This data would be a very 


valuable aid in efficient urban planning. 


NASA has been involved in the development one 
instrumentation to perform the variety of sensing required 
and demonstration of the data available to potential users. 
In addition, NASA has been asked by potential users to 
estimate total program cost and to recommend base locations 
for a special fleet of aircraft to support a program of 
remote sensing of urban centers. Assuming an initial cost 
for the program set up and instrument procurement and a 
fixed cost associated with each sensing mission, NASA has 
been seeking methods of analyzing the reduction of operating 
costs through efficient vehicle routing and effective base 


locations. 


The purpose of this paper is to develop an algorithm 
which will, given a base location and a set of urban areas 
to be sensed, provide an optimal or near optimal routing for 
the number of tours required to visit each of the areas. It 
1s assumed that the operating cost will be primarily a 
factor of distance traveled plus a fixed charge for each 
tour required to visit the set of urban areas. Using the 
fewest flights (subtours) to visit each area and optimally 
routing the aircraft to insure the minimum distance traveled 


Will yield the minimum operational cost. 


This problem was formulated as a multiple traveling 
salesman problem with the addditional constraint of aircraft 
range. The aircraft or sensing platform may only be capable 
of visiting a subset of the set of urban centers before it 
must return to its base (completing a subtour) to refuel. 
Thus the problem is to minimize the sum of the distances of 


the subtours. The algorithm presented by Raymond [ref 14] 





was modified to consider this range constraint. The 
modified algorithm was programmed in FORTRAN and will be 


discussed along with operational results and recommendations 


for additional uses. 





The traveling salesman problem (TSP) is one of the 
classic challenges in operations research (OR). The probiem 
is simply stated: a salesman starting from his home wishes 
to visit a set of cities once and only once and return hone. 
In what order should the cities be visited to minimize the 
total distance traveled? Mathematical formulation of the 
problem anda discussion of the constraints can be found in 
Wagner {ref 17] Chapter 13. 


OR literature has included applications of this model to 
a variety of practical problems. Industrial organizations 
have used the model to improve the routing of delivery 
vehicles, an important consideration in light of increased 
Oil prices and rising driver salaries. The "Chinese postman 
problem" is a Similiar problem applied to postal service 
introduced by the Chinese mathematician Kwan [ref 10]. 
Municipal governments have found the model valuable for 
routing garbage collection, street cleaning, meter reading 
and snow removal vehicles. Other applications include dairy 
and newspaper delivery. Further discussion of applications 
can be found in Bradley's "Survey of Deterministic Networks" 
{ref 3]. 


The key difference in the classic statement of the fTSP 
and the remote sensing problem is the addition of a range 
constraint. The salesman of the TSP may refuel during his 
trip without altering the optimal route. This is not 
feasible for the aircraft performing remote sensing. There 
exists a pcint during the trip when it must return to its 


base to refuel and then continue the mission. Therefore, 
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the problem statement becomes: an aircraft starting from its 
base must visit a set of urban areas once and only once and 
return to Ats base. The aircraft is subject to a range 
limitation for any given flight; that is, it may only be 
capable of visiting a subset of the set of cities before it 
Mitseeceturn to its base to refuel. Each time the aircraft 
returns to its base a subtour is completed and additional 
flights (subtours) are created as necessary to accomplish a 
misit to each city. The objective is tc minimize total 


distance traveled as the sum of the length of each subtour. 


A subtour, t, can be representec as a set of ordered 
city pairs, e.d. 
tem a (2 yi ) « « « (2 sla ry Ce aes Oa 
nie) 3 3 p-i non 1 


where is always the base. The mathematical formulation 


of the problem is to find the set of subtcurs to: 


s 
minimize > PE P| 
k=1 
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each city in exactly one subtour 
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where Ss = number of subtours 


ae = length of subtour k 

a = number of cities in subtour k 
d= distance from city i to city j 
1) 


cC = range of aircraft 
nc = number of cities to be visited 


Note that if s=1, eq. 2.2 replaces eq. 2.1 and the problem 
1s a standard TSP. The upper bound of eq. 2.4 represents 
the worse case Situation that each city must be visited by a 


seperate subtour. Eq. 2.3 represents the range constraint. 


The distance matrix (D. ) represents the euclidean 
13 


distance from city i to city j. During actual sensing 


Missions the aircraft must fly several passes over a city to 
completely image the area. This distance (loiter distance) 
depends on three factors: the size of the area, the percent 
Side overlap used for film images and the turning radius of 
the aircraft. Each city will have a unique loiter distance 
er that must be considered when the subtour length is 


calculated. When the problem is formulated the value ld is 
il 


added to each element in column i. Thus the new distance 


taeeixeets De = D + id for alli. The new matrix is 
1j 1j aI 
used to compute subtour lengths. 


This protlem appears similiar to two existing expansions 
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of the TSP; the multisalesman problem (MTSP) and the vehicle 
dispatcning problem. In the MTSP, m salesmen starting from 
a base city must visit n cities and return. Each city is 
visited by one and only one of the m salemen. The objective 
is to route each of the m salesmen such that the total 
distance traveled 1S minimized. Bellmore and Hong [ref 2] 
provide a detailed discussion and formulation of this 
problem. The objective of the vehicle dispatching problem is 
to obtain delivery routes from a depot to a set of demand 
points for a fleet of vehicles such that the total distance 
traveled by the entire fleet 1S minimized as with the a 
traveling salesmen. The problem is enhanced by finite 
vehicle capacities, route time constraints and varying 
quantity demanded at each demand point. A complete 
description cf this type application can be found in Golden 
{ref 12]. 


There are important differences between the remote 
sensing problem and the applications mentioned above. The 
MTSP has no upper bound on the number of cities any one of 
the m salesmen can visit. The aircraft range constraint 
requires consideration of a upper bound on the number of 
Cities that can be visited during any given flight. The 
vehicle dispatching problem is concerned with routing each 
vehicle in the fleet. ‘Although potential users will 
presumably have a "fleet" of aircraft to accomplish its 
sensing missions, not every aircraft will be scheduled to 
sense a subset of the areas to be visited. For example, 
given a fleet of three aircraft and the determination that 
two subtours are required to visit the set of cities, only 
two aircraft (or one aircraft making two flights) will be 
scheduled. This decision assumes the fixed cost charged 
each time an aircraft is launched is sufficiently high to 
make it more desireable to create the fewest number of 


flights. 
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The remote sensing problem, therefore, is not a_ simple 
restatement of an existing TSP. The addition of a range 
constraint requires special considerations. The solution 
technigue must carefully create only the humber of subtours 
required and route each subtour so that the total distance 
traveled is minimized. 
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| One of the most intriguing aspects of the TSP is, that 
although simply stated, it remains unsolved ina closed 
form. The computational difficulty with known solution 
tecnnigues is the exponential increasing solution time as 
the size of the problem increases. There are two types of 
solution algorithms: heuristic and exact. The latter 
yields the optimal tour. The difficulty with this type 
algorithna is the large computer storage and time 
requirements which are related to the large number of 
possible solutions. For example, given a symmetric TSP 
there are (n-1)!/2 possible solutions (n=number of nodes). A 
relatively small problem of ten nodes yields 181,440 
solutions. Fortunately, many applied or real world problems 
do not demand the optimal solution. Powerful heuristic 
algorithms have been developed which will yield near optimal 
solutions and provide a substantial savings in computer 
storage, time and money. An excellent survey of existing 


solutions is given by Bellmore and Nemhauser {ref 2]. 


The Bellmore and Nemhauser article classifies the 
algorithms by soiution generation method. There are three 
fundamentally different ways of generating solutions: 
tour-to-tour improvement, tour building, and subtour 
elimination. The remote sensing problem requires tour 
length feasibility be maintained as the tour-is formed. 
That is, if visiting the next city (node) will result in 
violating the range constraint, a new subtour is started. 
Consequently, a tour building solution was considered the 


best technique to adapt to this problen. 
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Both exact and heuristic algorithms were studied for 
applicability. A literature survey provided a number of 
both types of algorithms. An exact solution considered was 
presented by Christofides and Eilon {ref 5}. Their paper 
discusses three solution methods to the vehicle-dispatching 
problen. The branch and bound approach, a modification of 
the algorithm published by Little et. al. {ref 11], would 
reach an exact solution to this problem provided the number 
of subtours reguired was known a priori. If s (the number 
of subtours) was known, the base location would be replaced 
Byes artificial bases as explained in ref 5. The algorithm 
would progress uSing the branch and bound method to insert 
nodes into one of the tours. Before branching to a new node, 
it would be necesSary to check that the total distance 
accumulated on each subtour did not exceed the range. The 
key to this algorithm is knowing the number of subtours a 


Berori. 


As noted in the problem formulation section, there 
exists an upfer bound on the number of subtours. Therefore, 
it is feasible to select an s known, a priori to be 
sufficiently large, and apply this branch and bound 
algoritha. The solution would be recorded, the number of 
Subtours decreased by one and a new soluton computed. This 
procedure would continue until there are insufficient 
Subtours to complete a visit to each city and the last 
recorded solution would be the optimal. The success of this 
procedure would be very dependent upon the initial selection 
of s. A fpoor selection could involve many iterations and 
would require considerable user interaction and computer 
time. This would be particularly undesireable for problems 
that require a solution in a reasonably quick time. 
(Problems of this type are discussed in Section VI). 
Therefore, application of this algorithm was not considered 
profitable. 
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Limiting the study of heuristic algorithms to those 
using the tcur building generation technique, the algoritha 
presented by Raymond [ref 14] was selected to modify and 
apply to this problem. The approach of Raymonds's algoritaa 
is similiar to that of Karp and Thompson {[ref 9] and will 
obtain optigwal or near optimal solutions to the classical 
TSP. The algorithm successively inserts nodes into the tour 
by using a series of decision rules and provides for an 
exchange procedure to determine if the total distance 
traveled may be reduced by rearranging the order in which 
the nodes are visited. A detailed explanation of the 


algorithm steps is available in ref 14. 


The modified algorithm starts by creating a subtour fron 
the base to the most distant city. The same criteria as the 
original algorithm is used to select the next node for 


insertion. That is, increment calculations (I are made 
for each remaining node (k). This calculation is made by 


inserting the node k between each pair of nodes (p,q) in the 
subtour and computing the incremental increase in subtour 


length, I =d +d - d . Each node will have a mininun 
k pk kq Pq 
increment calculation (an) and a next best increment 
calculation a: If the number of cities in the tour is 
less than four the maxinua aan determines the node k to be 
inserted between the nodes that produced is Otherwise, 
the node k for which . = Le is a maximum is selected and 
inserted between the nodes that produced se a Before 


completing the insertion, the algorithm looks ahead to 


determine if the range constraint will be violated. If it 
1s not violated the node is inserted and the three Link 


exchange procedure described by Raymond is applied. fThis 
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procedure takes each node in turn out of its present 
subtour, reevaluates its insertion increments in that 
subtour and inserts the node between the nodes that now give 
the Minimum tour length increment. If the constraint is 
violated the node is not inserted and an additional subtour 


is formed. 


The nearest city to the base not yet in a subtour is 
used to form the next subtour. The criteria for inserting 
nodes for the multiple tour problem is changed. A node is 
inserted in the subtour that will create the least increase 
in total distance (sum of subtour lengths). A procedure was 
developed that exchanges nodes among subtours to determine 
if a savings in total distance traveled is possible by 
placing a node in a different subtour. Before any node 
insertions or exchanges among subtours is completed, the 
algorithm looks ahead to insure tour length feasibility is 
maintained. -The algorithm progresses iteratively until each 
city is visited. ) 
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IV. THE ALGORITHM 


The steps of the algorithn: 


1. With the base as one node select the node farthest fron 
base (assumed to be j*) and form a tour with links (3,7j*) 
and (j*,B). Set the number of subtours Ss equal to 1 and the 


number of links - equal to 2. 


2. Define an as the minimum tour-length increment that 
results from the insertion of node k in each existing link 


(p,q). Define nodes p* and q* as the nodes connected by the 
link (p*,q*) that produced Le In addition define ae as 


the increment for the node kx most nearly equal to ae Thus 
for node k, Ie 1s the best increment, ean the next best 
increment and (p*,q*) the link to break if inserting node k. 
Calculate ae and 7 for each node that is not currently 


included in the tour. 


a. ea Be < 4 choose the node a for which nan is a maximum. 
igi - > & choose the node mn for which aes = aa is a 
maximum. 


4. Insert the node m into the tour between p* and q* by 
adding links (p*,m) and (m,q*) and deleting link (p*,q*). 
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Calculate the new tour distance ae sta ar < range, complete 
the insertion and increase the number of links in the tour 
by one. If or > range, do not insert node m and go to step 


7 to create a new subtour. 


5. Take each node in turn out of the present tour, 
reevaluate its insertion increments and insert the node 
between the pair of nodes that now gives the minimum tour 
length increment. (In most cases the node is reinserted in 
the link created when the node was taken out of the tour, 


producing no change.) 


6. If all nodes have been added to the tour stop, otherwise 
go to step 2. 


7. Select the node m*, not already in a subtour, closest to 
the base and form a tour with links (B,m*) and (m*,B). Set 
the number of subtours s equal to st1 and the number of 


links 1 equal to 2. 
s 


8. Calcwate ate aS in step 2 for each node k remaining to 
be added to a subtour. Consider the subtours in the order 


in which they were created. For each subtour select the 


node that produced the mininun a. Before inserting the 
node two checks are made. First insure that the range 
constraint will not be violated. Second insure that the 


total distance traveled (sum of subtour lengths) will not be 
less if the selected node is inserted in a different 
Subtour. Only those nodes which pass these checks are 
inserted. Thus at any given iteration a node need not be 


inserted in each subtour. A given node may have a smaller 
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aed in one subtour, but if inserted, would force range 
constraint violation. In this situation place the node in 


the next best subtour. If nodes may not be added to any 
existing subtour without violating the range constraint go 


to step 7. 

9. Repeat step 5 for each subtour to which a node was added. 
10. For each subtour make the increment calculations using 
the nodes from the other subtours. Select the minimum I* 


and check if the total distance traveled may be reduced by 


adding the selected node to the given subtour and deleting 
it from its original subtour. If the distance is reduced 
make the change and repeat step 5. 


11. If all nedes are contained in a subtour stop, otherwise 
go to step 8. 
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The computer program was used to provide routes fron 
bases at Moffett Field, California, Kansas City and Wallops 
Station, Maryland to fifty selected urban centers. Figure 1 
is the grid system used to calculate the distances between 
cities. The cities were placed in one of three sections by 
the following criteria: the straight lines from Moffett 
Field to Kansas City and Wallops Station to Kansas City were 
bisected by a perpendicular, as shown in Figure 1, creating 
three sections. The cities would be imaged from the base in 
their section. Table 1 provides the coordinates and loiter 
distance of each city. The range of the aircraft was 
assumed to be 2900 statute miles. Tables 2-1, 2-2 and 2-3 


are the routing recommendations of the computer progran. 


Tables 3-1, 3-2 and 3-3 are the results of a feasible 
solution to the same problem computed at NASA's Ames 
Research Center. The constraints observed for determining 
this solution were: the aircraft would have a range of 2900 
statute miles, no more than three urban centers would be 
imaged during any subtour and flights from Moffett Field and 
Wallops Station would be held to a reasonable minimun. 


The NASA solution required 29 subtours anda total 
distance of 52,660 statute niles. The algorithm solution 
constructed 15 subtours ‘with a total distance of 36,505 
Statute miles traveled. A savings of 16,155 statute miles 
Cesults from the algorithm solution. This large savings is 
due, in part, to the differing criteria used to formulate 
the solutions. However, it does show the usefulness of the 
algorithm and the speed at which the program can compute a 


final solution. 
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TABLE 1 


CITY COORDINATES AND LOITER DISTANCE 


City 
Albany 
Allentown 
Atlanta 
Baltimore 
Birmingham 
Boston 
Buffalo 
Chicago 
Cincinnati 
Cleveland 
Columbus 
Dallas 
Dayton 
Denver 
Detroit 
Greensboro 
dart fond 
Houston 
Indianapolis 
Kansas City 
Los Angeles 
Louisville 
Memphis 
Miami 


Milwaukee 


X 
2230 
2190 
less, 
2190 
1820 
2350 
2030 
1660 
1830 
1930 
1900 
1350 
1820 

910 
1850 
2150 
2300 
1470 
1730 
1380 

220 
1780 
1560 
2300 
1620 
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s4 
1390 
1230 
700 
17150 
650 
14 70 
13 10 
1120 
1015 
1160 
1080 
500 
1080 
900 
1230 
910 
1350 
350 
1030 
880 
560 
910 
710 
370 


120. 


Loiter 
Distance 


190 
190 
230 
340 
340 
280 
240 
630 
350 
410 
230 
740 
250 
450 
B50 
290 
150 
670 
400 
380 
670 
190 
220 
370 
270 


Loiter distances provide by Ames Research Center. 





Git y 
Minneapolis 
Nashville 
New Orleans 
New York 
Norfolk 
Oklahoma City 
Omaha 
Philadelphia 
Phoenix 
Pittsburgh 
Portland 
Providence 
Rochester 
sacramento 
St Louis 
St Petersburg 
Salt Lake City 
San Antonio 
San Bernardino 
San Diego 
San Francisco 
Seattle 
Syracuse 
Toledo 
Washington DC 
Moffett Field 


Kansas City 


Wallops Station 


TABLE 1 


X 
1400 
1780 
1740 
2300 
2270 
1290 
1300 
2230 

560 
2030 

150 
2360 
2090 

150 
1590 
2140 

5 80 
1270 

310 

320 

100 

190 
2160 
1830 
2200 

120 
1380 
2270 
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(con! d) 
CITY COORDINATES AND LOITER DISTANCE 


4 
1260 
800 
410 
1300 
2020 
660 
1010 
1200 
4 70 
1150 
1250 
1370 
1330 
850 
370 
410 
910 
270 
560 
480 
780 
1420 
1350 
1150 
1100 
760 
880 
1080 


Loiter 
Distance 


380 
270 
290 
730 
230 
330 
230 
480 
250 
430 
450 
160 
150 
410 
500 
250 
200 
290 
190 
510 
530 
500 
310 
230 
350 
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FLIGHT MISSIONS FOR URBAN CENTER COVERAGE 
Bact:  MOEFETIT FIELD 


btour Cities Imaged Length 
Ly Sacramento - Seattle - Portland 2892 
| 2 Los Angeles - San Diego - San Bernadino 2184 

3 San Francisco - Salt Lake City - Phoenix 2696 


Total Distance = 7772 miles 


Time to Compute Solution= 1.27 seconds 
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AR BHo rr are 


FLIGHT MISSIONS FOR URBAN CENTER COVERAGE 
BASE: KANSAS CITY 


btour Cities Imaged 

4 Houston - San Antonio - Oklahoma City 

2 Dallas - Denver 

3 Omaha 
4 Indianapolis - Louisville - St. Louis - Kansas 
5 Nashville - Birmingham - New Orleans - Memphis 
6 Minneapolis - Milwaukee - Chicago 


Total Distance = 14092 


Time to Compute Solution = 6.97 seconds 


Pe 


City 


Length 
2878 
ZBo4 

581 
Zo 
2746 


2574 





DAB GE 3 


PLIGHT MISSIONS FOR URBAN CENTER COVERAGE 
BASE: WALLOPS STATTION 


ubtour Cities Imaged Length 
1 St Petersberg - Miami - Norfolk 2642 
2 Columbus - Cincinnati - Dayton - Detroit - Buffalo 2816 
3 Hartford - Providence - Boston - Albany - Philadelphia 2401 
4 New York - Syracuse - Rochester - Allentown - Baltimore 2607 
5 Washington, DC - Pittsburgh - Cleveland - Toledo 2451 
6 Greensboro - Atlanta 1724 


Total Distance = 14641 


Time to Compute Solution = 17.02 seconds 
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Subtour 


Ve, 


Tap LG. o= | 


FLIGHT MISSIONS FOR URBAN CENTER COVERAGE 


Base Ore slr FLELD 
Cities Imaged 

Salt Lake City 

Phoenix 

Denver - Omaha 

Houston 

San Antonio 

Dallas 

Oklahoma City - Kansas City 

New Orleans - Memphis 

St Louis - Minneapolis 

Chicago - Milwaukee 

Miami 

St Petersburg 

Birmingham - Atlanta 

Nashville - Louisville - Indianapolis 

Cincinnati - Dayton - Columbus 

Detroit - Toledo 

Cleveland - Pittsburgh 

Buffalo —- Rochester 


Syracuse - Albany 


Total Distance = 39520 miles 


29 


Length 


50 


2340 


1990 


2020 


13.90 


130 


1410 


1960 


1850 


1800 


2780 


2250 


2070 


ee) 


2100 


1880 


2380 


Zoe 


2710 





TAB EE) 3—2 


FLIGHT MISSIONS FOR URBAN CENTER COVERAGE 
BASE: KANSAS CITY 


Subtour Cities Imaged Length 
1 Seattle 1950 
2 Portland 1450 
3 San Francisco - Sacramento 1140 
4 Los Angeles - San Bernadino 160 
5 San Diego 1400 


Total Distance = 7450 miles 
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FABLE 3-3 


FLIGHT MISSIONS FOR URBAN CENTER COVERAGE 
BASE: WALLOPS STATION 


Subtour Cities Imaged Length 
1 Boston - Providence - Hartford 1430 
Zz New York 1120 
3 Philadelphia - Allentown 1080 
4 Washington, DC - Baltimore 970 
> Norfolk - Greensboro 1090 


Total Distance = 5690 miles 
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The purpose of this paper was to develop an algorithno 
that will provide optiimal or near optimal routing for a 
range constrained vehicle visiting a set of nodes. The 
algorithm may be used to solve similiar problems. Several 
considerations for use of this algorithm are discussed in 
this section. They are: under-flying satellite coverage, 


isolated disaster coverage and point reconaissance. 


Satellites are often used to image large land masses and 
bodies of water. These orbiting platforms perform the 
sensing misSion during periodic passes over the critical 
area. The product of a mission may sometimes be impaired by 
Cloud cover or may reveal the need for finer resolution 
imagery. The areaS requiring additional coverage are 
assumed dispersed and relatively small. These areas can be 
considered as the cities of the original problem and the 
algorithm used to determine, from a given base, the number 
of flights required, the optimal or near optimal routing and 


the total distance traveled. 


Disasters such as tornados or flooding rcivers may 
require government agencies to estimate damage and nake 
recommendations for appropriating the relief effort. Remote 
sensing can provide valuable data to assist in the decision 
process. Consider a flooding river; it is likely several 
areas along the flood path need immediate attention. The 
goal is to provide the data for evaluation as rapidly as 
possible. The algorithm can provide routing from a base to 
the areas of concern such that the sensing mission is 


completed in the shortest possible time. It is assumed that 


32 








least distance routing implies the least time will be 


required to accomplish the mission. 


Several other applications are similiar to the disaster 
coverage and involve routing reconnaissance aircraft. lo a 
example, when fighting large forest fires, “hot spots" 
develop that need periodic monitoring. Also, in a large 
theater of operations, the military commander often requires 
timely reconnaissance of key objective areas. In both these 
Situations the goal is to provide the information to the 
user aS rapidly as possible. Therefore, the algorithm can 
be used to route the aircraft to the critical areas and 


provide the information as quickly as possible. 
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Wie ooONc LUSLONS 


The algoritha and computer program presented can 
successfully construct an optimal or near optimal solution 
to the range constrained vehicle routing problem formulated 


in this paper. 


Several aspects of remote sensing, considered beyond the 
scope of this paper, provide a basis for further study in 
this area. The angle of the sun can have an appreciable 
effect on the resolution of the imagery. Sun-angle is a 
function of the time of day and could be considered in the 
computer program by adding a time dimension. Using aircraft 
speed, it would be possible to calculate a time period 
during which the aircraft would image a given city. This 
would provide the user the ability to insure that the 


desired sun-angle occurs during the imaging. 


This paper assumes each aircraft has sufficient sensor 
Capacity to image any number of cities in a subtour. It is 
likely there may be a finite capacity for certain sensors. 
For example, assume only a fixed amount of film can be 
placed on the aircraft. The film requied to image all the 
cities ina given subtour may exceed this amount. To 
consider this problem, the aircraft would start with a 
finite amount of film which would be reduced by the amount 
required to image a given city as that city is visited. 
Thus sensor capacity could be added as a constraint and 


considered in the same manner as the range constraint. 


The problem of imaging areas as widely dispersed as the 


fifty urban centers listed in Table 1 requires multiple 
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bases. The example problem assigned cities to bases by 
section. This approach is reasonable, however, it would be 
preferable to permit the computer program to make this 
asSignment. One consideration is to construct an artificial 
base from which all flights originate and sufficient dunny 
bases at the same location as the actual bases to permit 
multiple subtours through that base. The distance matrix 
must be constructed such that each subtour starts from the 
artificial base to one of the actual bases and is’ routed 
back through the actual base to the artificial base. 
Formulation of a fifty city problem with artificial and 


dummy bases could require considerable computer storage. 
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APPENDIX A 


COMPUTER PROGRAM 


WR WK RE HE eA CE GH HE RC ER EE BC EEK BK AE Oe KE EK AR AE OK i ee HE OK I OK AE OK OK oe eK 


HHH HHHHHHHKHHHHHHHHHHHHHHHKHHHHHHHKHHHHHKHHHHHHHHEHH EH 


VARIABLE DEFINITIONS 


X(I)=X COORDINATE OF THE I-TH CITY 

Y(I)=Y CQORDINATE OF THE I-+TH CITY 

Z(1)=LOITER DISTANCE OF THE I-TH CITY 
D(I,J)=DISTANCE MATRIX 

RANGE=AIRCRAFT RANGE 

NTGURS=NUMBER OF SUBTOURS CREATED 
LINKS(I)=NUMBER OF LINKS IN SUBTOUR I 
INNODE(IT,J)J=VECTOR OF CITIES IN SUBTOUR I 
NONODE(TI=VECTOR OF CITIES TO BE VISITED 
LENGTH(IJ=LENGTH OF SUBTOUR I 

TOTDIS=TOTAL DISTANCE TRAVELED 

LH(I,J)=LEFT HAND ENDPOINT OF SUBTOUR I, LINK J 
RH(I,J)=RIGHT HAND ENDPOINT OF SUBTOUR I, LINK J 
NC=NUMBER OF CITIES IN THE PROBLEM 

S(Ty,JI=IN SUBTIUR IT THE SUCCESSOR OF CITY J 
PRED(I,JIJ=IN SUBTOUR [I THE PREDECESSOR OF CITY J 


STAR (I yJ) > 

NNSTR( tJ) $ 

LOSHRICLT. J) $ 

SEeSuR (1 od) s STORAGE MATRICES FOR ALGORITHM 
SNNSTR (I,J) S CALCULATIONS 

SLLSTR(I+ J) $ 

THETA( I,J) > 


HHH HHHHHHHHKHHHHHHHHKHHHHHHHRHHHHHHHHHHHHHHHKH HH 


SESE LSCCOC LOSS STS SSCOOCOS SSCS TCSCSCoOSL SSS CCS SCL CL SS Ses 2. 


MAIN 
COMMON D( 25,25) 4S(25925) 9 INNODE (25,25) »NONODE (25), 
1LRH(25 925) o LHC 255 2549L INKS ( 25) THATA(25 425) gNTOURSe KKy 
LIEND NC y TOTDIS, LENGTH (25) sRANGEs ISTOP, PRED (25,25) 
DIMENSION X(25)+Y(25) 9 STAR( 254525) @NNSTR( 25925) 4 
LLESTR(25525) » SECSTRI25 925) 9 SNNSTR( 25925) g SLLSTR( 25925) 
1 sCRIT(25) «OLD 2H(25 425) pOLDLH(25,25) 9225) 
INTEGER S»RHy OLDRH,OLDLH» EXNODE » PRED 
REAMMEENG TH 
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